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Spontaneous formation of geometric patterns 
is a fascinating, ubiquitous process that provides 
fundamental insights into the roles of symmetry 
breaking, anisotropy and nonlinear interactions 
in emergent phenomena mm- Here we report 
dynamic, highly ordered topographic patterns on 
the surface of diamond that span multiple length 
scales and have a symmetry controlled by the 
chemical species of a precursor gas used in elec¬ 
tron beam induced etching (EBIE). This behavior 
reveals an underlying etch rate anisotropy and an 
electron energy transfer pathway that has been 
overlooked by existing EBIE theory. We present 
an etch rate kinetics model that fully explains our 
results and is universally applicable to EBIE. Our 
findings can be exploited for controlled wetting, 
optical structuring and other emerging applica¬ 
tions that require nano and micro-scale surface 
texturing. 

Electron beam induced etching (EBIE) [H |5] is a high 
resolution, single step, direct-write nanofabrication tech¬ 
nique in which a precursor gas and an electron beam 
are used to realize etching. To date, EBIE has been 
used to machine a wide range of materials using etch 
precursor such as oxygen, water, ammonia, nitrogen tri¬ 
fluoride, xenon difluoride and chlorine. Key advantages 
of EBIE include cite-speciflcity and the ability to etch 
materials such as diamond which are resistant to conven¬ 
tional chemical etch processes. Consequently, EBIE has 
recently been used to realize practical device components 
for use in photonics [6], plasmonics [7] and nanofluidics 

i. 

In this work, we report an emergent pattern formation 
phenomenon caused by a chemical etch rate anisotropy 
in EBIE of single crystal diamond. The results reveal a 
shortcoming in existing, established EBIE theory which 
does not adequately explain the observed etch kinet¬ 
ics. We therefore propose a fundamental modification, 
whereby the critical role of energetic electrons is to trans¬ 
fer energy to and break bonds between surface atoms 
of the solid rather than to surface-adsorbed precursor 
molecules. The new EBIE model is confirmed exper¬ 
imentally, explains the observed patterns, and resolves 
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long standing problems that have been identified in the 
EBIE literature. 

Eigure [^a) is a schematic illustration of EBIE per¬ 
formed using H 2 O precursor gas. Eigure [^b) shows im¬ 
ages of topographic patterns that form on the surface of 
single crystal (001) oriented diamond during H 2 O EBIE 
performed at room temperature. A movie showing the 
pattern formation and evolution dynamics is provided as 
Supplementary Video ffl. Such patterns have not been 
reported previously because their formation is inhibited 
by small amounts of residual hydrocarbon contaminants 
present on the etched surface. Residual contaminants are 
very common in electron microscopes, alter the surface 
termination during EBIE and give rise to a competing 
process of electron beam induced deposition of carbon 

0 . 

Etching initiates at scratches and other surface defects, 
which expand laterally during EBIE, evolving into highly 
symmetric rhombohedra such as the large pit seen in the 
top left corner of Eigure Bb). Similarly, the topogra¬ 
phy that is normally associated with surface roughness 
caused by EBIE rapidly evolves into step edges with 
{110} sidewalls which propagate laterally until they reach 
the edge of the area scanned by the electron beam. The 
{111} family of planes is absent from the resulting surface 
topography, the step sidewalls are comprised of {110} 
planes, and 90° step corners are formed at the intercepts 
of the {110} planes. Corner formation requires the (110), 
(IlO), (110) and (110) planes to etch slower than the 
(100), (010), (010) and (lOO) planes. Erom these obser¬ 
vations, we conclude that H20-mediated EBIE removes 
material from the {100} and {111} planes faster than 
from the {110} family of planes. 

In order to prove conclusively that the proposed 
anisotropy yields the four-fold symmetry observed in Eig¬ 
ure Bb), we used the 3D implementation of the level set 
method (LSM) pT] . LSM is a robust technique for evolv¬ 
ing implicit surfaces under anisotropic velocity fields. It 
can be used to calculate surface shapes produced by etch 
rate anisotropies defined by differences between the etch 
rates of specific crystal planes. The simulations detailed 
in the Supporting Information reveal that the calculated 
surfaces match experiment only if the {110} planes that 
are oriented at 90° with respect to the electron beam axis 
etch slower than all other planes. Supplementary Video 
ff2 and Eigure [^d) show the resulting surface features 
(rhombohedra). 

Validity of the simulation was confirmed by applying 



2 


the same etch rate anisotropy rule set to ( 111 ) oriented 
diamond. The simulation and H 2 O EBIE both produced 
the trigons shown in Eigurej^c) and (e). The LSM simu¬ 
lations therefore support our conclusion that the geome¬ 
tries of patterns observed during H 2 O EBIE of (001) and 
( 111 ) oriented diamond is governed primarily by slow 
etching of specific {110} planes. This anisotropy can 
not be explained by conventional, established EBIE the¬ 
ory which is based on the mechanistic framework shown 
schematically in Eigurej^a), where the key role of ener¬ 
getic electrons is to dissociate surface-adsorbed precursor 
molecules. In the case of H 2 O EBIE of diamond, a pos¬ 
sible pathway in this framework is the following Il[5]: 
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where the subscripts [v], [s], [p] and [c] signify vapor 
phase, solid phase, physisorbed and chemisorbed species, 
respectively. represents the energy barrier for dissoci¬ 
ation of H 2 O, and E 2 is the binding energy of the reaction 
product. According to the standard EBIE model, 5i and 
E 2 are overcome by a transfer of kinetic energy from the 
electrons that drive EBIE, and thermal energy of the sub¬ 
strate (/cT), respectively. This model has been used to 
explain a wide range of experiments such as dependencies 
of etch rates on time, beam current density and pressure 
of the precursor gas 00 [TTH 2 T] . However, the model 
can not explain the etch rate anisotropy seen in Eigure[^ 
unless different crystal planes give rise to significant vari¬ 
ations in the electron dissociation cross-section of H 2 O 
adsorbates, the secondary electron emission yield, or the 
local coverage of precursor molecule adsorbates. None of 
these are plausible since the precursor is H 2 O, the pat¬ 
terns form at room temperature, and the slowest etching 
planes are not consistently dark in secondary electron 
images. 

To resolve the above issues, we propose a new mecha¬ 
nism, in which electrons provide the energy E 2 in Reac¬ 
tion]^ as shown in Eigurej^b). That is, the critical role 
of electrons is not to dissociate the physisorbed precursor 
molecules, but to break bonds that bind surface atoms 
to the substrate and thus enable the desorption of the 
final reaction products. During etching, Reactionj^is ex¬ 
pected to proceed spontaneously since active surface sites 
are generated continuously and the precursor molecules 
will dissociate on unterminated sites. In this framework, 
the etch rate anisotropy needed to produce the patterns 
seen in Eigure ^is expected since ^2 (i.e. the C-C bond 
strength and the corresponding cross-section for scission 
by electrons) varies with the crystal plane. To confirm 
the proposed mechanism we performed an experiment 
based on the fact [ 22 ] that the C-C bond strengths are 
modified by hydrogen which reconstructs and stabilizes 
the {111} surface. We therefore performed H 2 O EBIE of 


( 001 ) and ( 111 ) oriented diamond in the presence of NH 3 
gas, where the role of the NH 3 is to supply an excess of 
hydrogen radicals to terminate the ( 111 ) planes. Eigure 

shows that the corresponding surface patterns consist 
of inverted pyramids and trigons, respectively, and that 
these geometries are indeed expected from LSM simula¬ 
tions in which the { 111 } planes are the slowest etching 
planes. 

We performed one more experiment to further test the 
proposed EBIE mechanism. A consequence of the con¬ 
ventional EBIE model is that the EBIE rate is directly 
proportional to the concentration of physisorbed precur¬ 
sor molecules 00 . Hence, the etch rate of diamond is 
expected to be negligible at a temperature of ^ 400 K, 
as is shown in Eigure]^ (solid curves), irrespective of the 
electron flux used to perform EBIE. However, we ob¬ 
serve significant etch rates at temperatures as high as 
600 K. This result can not be explained by the conven¬ 
tional model, but is consistent with the new model in 
which the EBIE rate is proportional to the concentration 
of chemisorbed oxygen. The observed temperature de¬ 
pendence therefore serves as direct evidence for the new 
EBIE model. 

We note that the new EBIE model (Eigure [^b)) is 
consistent with all results in the literature that were ex¬ 
plained successfully by the conventional EBIE model. 
Eirst, both models predict the same dependence of EBIE 
rate on electron beam energy since the secondary elec¬ 
tron yield has the same dependence on beam energy as 
the amount of energy that is deposited by the beam into 
the surface atoms of the substrate. Second, both mod¬ 
els predict the existence of reaction rate limited and mass 
transport limited etching regimes (as is shown in the Sup¬ 
porting Information), which makes both models consis¬ 
tent with a large amount of experimental data available 
in the EBIE literature. However, the new model is unique 
in being consistent with reports of UV laser induced etch¬ 
ing of diamond, that is believed to proceed through a two 
photon C-C bond scission mechanism [3]. Eurthermore, 
the new model provides a satisfactory explanation for the 
fact that single crystal diamond can be etched by EBIE 
in the first place. The energy barrier of Reactionj^in dia¬ 
mond is known to be significant [23] and therefore etching 
observed at room temperature, or any temperature below 
the onset of defect generation and graphitization cannot 
be accounted for in the standard model. 

Einally, we note that the topographic patterns cannot 
be explained by an anisotropic sub-surface damage gen¬ 
eration mechanism analogous to the graphitization path¬ 
ways encountered in conventional dry and wet diamond 
etch processes [24[ |25]. Eirst, the etch rate anisotropy 
was modified significantly by the presence of NH 3 gas, 
which should not change the sub-surface damage gen¬ 
eration rate. Second, prior studies of EBIE of single 
crystal diamond have failed to produce any evidence of 
damage by photoluminescence and Raman spectroscopy 
[SlEollMlEZ]. Third, the generation rate of damage pro¬ 
duced by a 5 keV electron beam scales with the local 
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energy density deposited into the substrate throughout 
the electron interaction volume [28] . The damage gener¬ 
ation rate is therefore isotropic, except for special cases 
where the electron beam axis is parallel to a channel¬ 
ing axis, which should produce a strong dependence of 
the patterns on sample tilt, which was not observed in 
our experiments. We therefore conclude that sub-surface 
damage generation does not play a role in the observed 
etching and pattern formation behavior. 

To summarize, we showed several dynamic pattern for¬ 
mations on the surface of single crystal diamond. We 
proposed an amended model for the EBIE process that 
is based on interactions of electrons with the substrate 
rather than the precursor molecule adsorbates. The new 
model explains our results and rectifies long standing is¬ 
sues in the literature of EBIE. Our results can be lever¬ 
aged to engineer surface patterns controlled by electron 


beam irradiation conditions. 
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FIG. 1. Topographic patterns formed during H 2 O mediated electron beam induced etching of single crystal 
diamond, (a) Schematic illustration of H 2 O EBIE. (b) Expanding rhombohedra formed on the surface of (001) oriented 
diamond, and (c) trigons on the surface of (111) diamond. (d,e) Corresponding simulated rhombohedra and trigons (colored 
by the relative local etch rate) that are expected if the {110} planes are the slowest etching planes. 
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EIG. 2. Simplified schematics of the mechanisms behind the standard and new EBIE models, (a) Standard model: 
(i) Electrons dissociate physisorbed precursor molecules generating radicals that (ii) react with surface atoms and (iii) desorb 
spontaneously at the temperature used to perform EBIE. (b) New model: (i) Gas phase and physisorbed precursor molecules 
decompose spontaneously, and (ii) electrons break bonds in the crystal lattice causing (iii) the desorption of etch reaction 
products. 
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FIG. 3. Topographic patterns formed during electron 
beam induced etching of single crystal diamond in the 
presence of NHg. (a) Expanding inverting pyramids formed 
on the surface of (001) oriented diamond, and (b) trigons on 
the surface of (111) diamond. (c,d) Corresponding simulated 
pyramids and trigons (colored by the relative local etch rate) 
that are expected if the {111} planes are the slowest etching 
planes. 


10 


^ 10 ' 


g 10' 


10 


10 ^ 



250 300 350 400 450 500 550 600 650 


Temperature (K) 


FIG. 4. Temperature dependence of the rate of EBIE. 

Rate of EBIE measured experimentally (points), and calcu¬ 
lated using the established model of EBIE detailed in the 
Supporting Information (lines) using a wide range of electron 
fluxes. The etch rates are normalized to the rate of EBIE at 
room temperature. 








